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A search for the rare decays B® — > yT and B — > fi + fi~ is performed using data collected in 
2011 and 2012 with the LHCb experiment at the Large Hadron Collider. The data samples comprise 
1.1 fb _1 of proton-proton collisions at yfs = 8 TeV and 1.0 fb _1 at yfs = 7 TeV. We observe an excess 
of B° — > candidates with respect to the background expectation. The probability that the 

background could produce such an excess or larger is 5.3 x 10~ 4 corresponding to a signal significance 
of 3.5 standard deviations. A maximum-likelihood fit gives a branching fraction of B(Bg — > fj, + fj,~) 
= (3.2 ti'2) x 10~ 9 , where the statistical uncertainty is 95 % of the total uncertainty. This result is in 
agreement with the Standard Model expectation. The observed number of B° — > fi + /i~ candidates is 
consistent with the background expectation, giving an upper limit of B(B° — > fi + < 9.4 x 10 -10 
at 95 % confidence level. 
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The rare decays B° s — > and B° —> ^t + /i~ are 

highly suppressed in the Standard Model (SM). Pre- 
cise predictions of their branching fractions, B(B° S — > 
fi+fj,-) = (3.23 ± 0.27) x 1CT 9 and B{B° -> = 
(1.07 ± 0.10) x 10~ 10 H] make these modes powerful 
probes in the search for deviations from the SM, espe- 
cially in models with a non-standard Higgs sector. Tak- 
ing the measured finite width difference of the B® sys- 
tem [2] into account [3], the time integrated branching 
fraction of £?° — » that should be compared to the 

experimental value is (3.54 ± 0.30) x 10~ 9 . 

Previous searches [U-[S] already constrain possible 
deviations from the SM predictions. The lowest 
published limits are B(B® — > < 4.5 x 10~ 9 and 

B(B° -> n + n~) < 1.0 x 10~ 9 at 95% confidence level 
(CL) from the LHCb collaboration using 1.0 ftT 1 of data 
collected in pp collisions in 2011 at y/s = 7 TeV [8]. This 
Letter reports an update of this search with 1.1 fb _1 of 
data recorded in 2012 at yfs = 8 TeV. 

The analysis of 2012 data is similar to that described in 
Ref. [5] with two main improvements: the use of particle 
identification to select B9-. —> h + h'~ (with = K,tt) 
decays used to calibrate the geometrical and kinematic 
variables, and a refined estimate of the exclusive back- 
grounds. To avoid potential bias, the events in the signal 
region were not examined until all the analysis choices 
were finalized. The updated estimate of the exclusive 
backgrounds is also applied to the 2011 data [8] and the 
results re-evaluated. The results obtained with the com- 
bined 2011 and 2012 datasets supersede those of Ref. [8]. 

The LHCb detector is a single-arm forward spectrom- 
eter covering the pseudorapidity range 2 < -q < 5, and is 
described in detail in Ref. [pj . The simulated events used 
in this analysis are produced using the software described 
in Refs. [MB]. 

Candidate B^ —> n + events are required to be se- 
lected by a hardware and a subsequent software trigger. 
The candidates are predominantly selected by single and 
dimuon triggers [T7] and, to a smaller extent, by a generic 
b- hadron trigger [TS]. Candidate events in the B + — » 
J/ipK + control channel, with J/tp — > (inclusion 
of charged conjugated processes is implied throughout 



this Letter), are selected in a very similar way, the only 
difference being a different dimuon mass requirement in 
the final software trigger. The 5°^ — > h + h'~ decays are 
predominantly selected by a hardware trigger based on 
the calorimeter transverse energy and subsequently by a 
generic 6-hadron software trigger. 

The -B° s ) — > n + p~ candidates are selected by requiring 
two high quality muon candidates [TH] displaced with re- 
spect to any pp interaction vertex (primary vertex, PV) , 
and forming a secondary vertex (SV) with a x 2 per degree 
of freedom smaller than 9 and separated from the PV in 
the downstream direction by a flight distance significance 
greater than 15. Only candidates with an impact param- 
eter x 2 , IPx 2 (defined as the difference between the \ 2 of 
the PV formed with and without the considered tracks) 
less than 25 are considered. When more than one PV 
is reconstructed, that giving the smallest IPx 2 for the 
B candidate is chosen. Tracks from selected candidates 
are required to have transverse momentum satisfying 
0.25 < p T < 40GeV/c and p < 500GeV/c. Only B can- 
didates with decay times smaller than 9t(-B 9 ) [2D] and 
with invariant mass in the range [4900, 6000] MeV/c 2 are 
kept. 

Dimuon candidates from elastic diphoton production 
are heavily suppressed by requiring p^(B) > 0.5GeV/c. 
The surviving background comprises mainly random 
combinations of muons from semileptonic decays of two 
different b hadrons (bb — > /i + /i~AT, where X is any other 
set of particles). 

Two channels, B + — > J/tpK + and B° — > K + tt~, serve 
as normalization modes. The first mode has trigger and 
muon identification efficiencies similar to those of the sig- 
nal, but a different number of tracks in the final state. 
The second mode has a similar topology, but is triggered 
differently. The selection of these channels is as close 
as possible to that of the signal to reduce the impact of 
potential systematic uncertainties. 

The B° — > K + ir~ selection is the same as for B? , — > 
signal except for muon identification. The two 
tracks are nevertheless required to be within the muon 
detector acceptance. 

The J/ij) -> decay in the B+ -> J/ipK + 



normalization channel is also selected similarly to the 
^(s) signals, except for the requirements on the 

IPX 2 and mass. Kaon candidates are required to have 
IPX 2 > 25. 

A two-stage multivariate selection, based on boosted 
decision trees [2T1 122] is applied to the B? s — s- fi + n~ 
candidates. A cut on the first multivariate discriminant, 
unchanged from Ref . [5] , removes 80 % of the background 
while retaining 92 % of signal. The efficiencies of this cut 
for the signal and the normalization samples are equal 
within 0.2 % as determined from simulation. 

The output of the second multivariate discriminant, 
called BDT, and the dimuon invariant mass are used to 
classify the selected candidates. The nine variables enter- 
ing the BDT are the B candidate IP, the minimum IPx 2 
of the two muons with respect to any PV, the sum of the 
degrees of isolation of the muons (the number of good 
two-track vertices a muon can make with other tracks in 
the event), the B candidate decay time, py, and isola- 
tion [23], the distance of closest approach between the 
two muons, the minimum of the muons, and the co- 
sine of the angle between the muon momentum in the 
dimuon rest frame and the vector perpendicular to both 
the B candidate momentum and the beam axis. 

The BDT discriminant is trained using simulated sam- 
ples consisting of B® —> for signal and bb —> 
for background. The BDT response is defined 
such that it is approximately uniformly distributed be- 
tween zero and one for signal events and peaks at zero 
for the background. The BDT response is independent 
of the invariant mass for signal inside the search window. 
The probability for a BP % — > event to have a given 
BDT value is obtained from data using B° —> K + tt~, 
7r + 7r~ and B® — > ir + K~, K + K~ exclusive decays se- 
lected as the signal events and triggered independently 
of the tracks from BP n candidates. 

The invariant mass lineshape of the signal events is 
described by a Crystal Ball function [24]. The peak 
values for the B® and B° mesons, m B o and m B o, are 
obtained from the B° -> K + K~ and B° -> K + ir~, 
B° — > n + ir~ samples. The resolutions are determined 
by combining the results obtained with a power-law in- 
terpolation between the measured resolutions of char- 
monium and bottomonium resonances decaying into two 
muons with those obtained with a fit of the mass distri- 
butions of B° -> K+7T-, B° -> 7T+7T- and B° -> K + K~ 
samples. The results are a B o — 25.0 ± 0.4MeV/c 2 and 
(j B a = 24.6 ± 0.4 MeV/c 2 , respectively The transition 
point of the radiative tail is obtained from simulated 
Bj? — > events smeared to reproduce the mass reso- 

lution measured in data. 

The B° — > n + jJT and B° — > fi + yT yields are translated 



into branching fractions using 

7-,n 4- —\ ^norm ^norm Jnorm at 

B(B° {S) -> M + M ) = — x N B o +fl - 

= "b™-> m + m - x ^Vb^ (!) 

where B norm represents the branching fraction, iV norm the 
number of signal events in the normalization channel ob- 
tained from a fit to the invariant mass distribution, and 
N B o ) _ S . /J + A1 - is the number of observed signal events. 

The factors fd( s ) and / norm indicate the probabilities 
that a b quark fragments into a B^ meson and into the 
hadron involved in the given normalization mode, respec- 
tively. We assume fa = f u and use f s /fd = 0.256 ±0.020 
measured in pp collision data at y/s = 7 TeV [25 . This 
value is in agreement within 1.5 a with that found at 
^fs = 8 TeV by comparing the ratios of the yields of 
B° — > J/ip<fi an d B + — > J/tpK + decays. The measured 
dependence of f s /fd on Pt{B) [2S] is found to be negli- 
gible for this analysis. 

The efficiency e s i g ( n0 rm) for the signal (normalization 
channel) is the product of the reconstruction efficiency of 
the final state particles including the geometric detector 
acceptance, the selection efficiency and the trigger effi- 
ciency. The ratio of acceptance, reconstruction and selec- 
tion efficiencies is computed using simulation. Potential 
differences between data and simulation are accounted 
for as systematic uncertainties. Reweighting techniques 
are used for all the distributions in the simulation that 
do not match those from data. The trigger efficiency 
is evaluated with data-driven techniques 26J. The ob- 
served numbers of B + — > J/ipK + and B° — > K + ir~ 
candidates in the 2012 dataset are 424 200 ± 1500 and 
14 600 ± 1100, respectively. The two normalization fac- 
tors a 1 ^" 1 j- are in agreement within the uncertain- 

ties, and their weighted average, taking correlations into 
account, gives aso_ +Al + M - = (2.52 ± 0.23) x 10~ 10 and 
auo_^ +/1 - = (6.45 S ±0.30) x 10~ n . 

In total, 24044 muon pairs with invariant mass be- 
tween 4900 and 6000MeV/c 2 pass the trigger and selec- 
tion requirements. Given the measured normalization 
factors and assuming the SM branching fractions, the 
data sample is expected to contain about 14.1 B° — > 
y + pT and 1.7 B° — > fi + decays. 

The BDT range is divided into eight bins with bound- 
aries [0.0, 0.25, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0]. For the 2012 
dataset, only one bin is considered in the BDT range 
0.8-1.0 due to the lack of events in the mass sidebands 
for BDT > 0.9. The signal regions are defined by 
m B o ±60MeV/c 2 . 

The expected number of combinatorial background 
events is determined by interpolating from the invariant 
mass sideband regions defined as [4900 MeV/c 2 , m B o — 
60 MeV/c 2 ] and [m B o + 60 MeV/c 2 , 6000 MeV/c 2 ] . The 
low-mass sideband and the B° and B° signal regions 
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are potentially polluted by exclusive backgrounds with 
or without misidentification of the muon candidates. 

The first category includes B° —> ir~ n + v^, B9s —> 
h + h'~, B° -t K^^+Vfj, and A£ -t ppTV^ decays. The 
B° —> n~ n + v^ and B9*. — > h + h'~ branching fractions 
are taken from Ref . [2D] . The theoretical estimates of the 
A" — > pii~v^ and B® — > K^fi+v^ branching fractions 
are taken from Refs. [27] and [28], respectively. The 
mass and BDT distributions of these modes are evalu- 
ated from simulated samples where the K — > fj,, ir — > fi 
and p — > fi misidentification probabilities as a function of 
momentum and transverse momentum are those deter- 
mined from D* + -» D a TT + 1 D° — > K~tt + and A — > pir~ 
data samples. We use the A 9 fragmentation fraction / A o 
measured by LHCb [2DJ and account for its pt depen- 
dence. 

The second category includes B+ — > J 
B® —> ^+^i _ 7 and B ^ — > 7r°( + '^ + ^~ decays, evaluated 
assuming branching fraction values from Refs. [3D], [3T] 
and [32) . respectively. Apart from B9^ — > h + h'~, all 
background modes are normalized relative to the B + — > 
J/ipK+ decay. The B° %~ n+v^, B? s) ->■ h + h'~ and 
— >. 7r°(+) fj,+ n~ decays are the dominant exclusive 
modes in the range BDT > 0.8, which accounts for 70 % 
of the sensitivity. 

In the full BDT range, 8.6 ± 0.7 doubly misidentified 
B9 , — > h + h'~ decays are expected in the full mass in- 
terval, 4.1±J;| in the 5° and 0.76±g;f| in the B° s sig- 
nal region. The expected yields for B° —> ir~ and 
B o(+) _^ 7r o(+) /i + Af - are 41 .i ±0 .4 and 11.9±3.5, respec- 
tively, in the full mass and BDT ranges. The contribu- 
tions of these two backgrounds above mgo — 60 MeV/c 2 
are negligible. The fractions of these backgrounds with 
BDT > 0.8, in the full mass range, are (19.0 ± 1.4)%, 
(11.1 ± 0.5)%, and (12.2 ± 0.3)% for B° s) ->■ h+h'~, 

B° —> 7r~/i + ^ t and B°( +s> — > n 01 --^ n + fi~ decays, respec- 
tively. 

A simultaneous unbinned maximum-likelihood fit to 
the mass projections in the BDT bins is performed on 
the mass sidebands to determine the number of expected 
combinatorial background events in the B° and B® sig- 
nal regions used in the derivation of the branching frac- 
tion limit. In this fit the parameters that describe the 
mass distributions of the exclusive backgrounds, their 
fractional yields in each BDT bin and their overall yields 
are limited by Gaussian constraints according to their 
expected values and uncertainties. The combinatorial 
background is parameterized with an exponential func- 
tion with slope and normalization allowed to vary. The 
systematic uncertainty on the estimated number of com- 
binatorial background events in the signal regions is de- 
termined by fluctuating the number of events observed 
in the sidebands according to a Poisson distribution, and 
by varying the exponential slope according to its uncer- 



tainty. The same fit is then performed on the full mass 
range to determine the B° — ¥ n + and B a — > \i + 
branching fractions, which are free parameters of the fit. 
The B a s — > and B Q — > fractional yields in 

BDT bins are constrained to the BDT fractions cali- 
brated with the B? s j — > h + h'~ sample. The parameters 
of the Crystal Ball functions that describe the mass line- 
shapes and the normalization factors are restricted by 
Gaussian constraints according to their expected values 
and uncertainties. 

The compatibility of the observed distribution of 
events with that expected for a given branching frac- 
tion hypothesis is computed using the CL S method [33j . 
The method provides CL s+ b, a measure of the com- 
patibility of the observed distribution with the signal 
plus background hypothesis, CLb, a measure of the 
compatibility with the background-only hypothesis, and 
CL S = CL s+ b/CLb- 

The invariant mass signal regions are divided into nine 
bins with boundaries m B o ± 18, 30, 36, 48, 60 MeV/c 2 . In 
each bin of the two-dimensional space formed by the 
dimuon mass and the BDT output we count the num- 
ber of observed candidates, and compute the expected 
number of signal and background events. 

The comparison of the distributions of observed events 
and expected background events in the 2012 dataset re- 
sults in p- values (1 — CL b ) of 9xl0~ 4 for the B" -> (i + fi~ 
and 0.16 for the B° — s- decay, computed at the 

branching fraction values corresponding to CL s+ b = 
0.5. We observe an excess of B® — > n + [iT candidates 
with respect to background expectation with a signifi- 
cance of 3.3 standard deviations. The simultaneous un- 
binned maximum-likelihood fit gives £>(£> 9 — > = 
(5.lt 2 ^(stat)±°^(syst)) x 10" 9 . The statistical uncer- 
tainty reflects the interval corresponding to a change of 
0.5 with respect to the maximum of the likelihood af- 
ter fixing all the fit parameters to their expected val- 
ues except the B° — > and B° —> branching 
fractions and the slope and normalization of the combi- 
natorial background. The systematic uncertainty is ob- 
tained by subtracting in quadrature the statistical un- 
certainty from the total uncertainty obtained from the 
likelihood with all nuisance parameters left to vary ac- 
cording to their uncertainties. An additional systematic 
uncertainty of 0.16x 10 -9 reflects the impact on the result 
of the change in the parameterization of the combinato- 
rial background from a single to a double exponential, 
and is added in quadrature. 

The expected and measured limits on the B° — > 
branching fraction at 90 % and 95 % CL are shown in 
Table [T] The expected limits are computed allowing for 
the presence of B?^ — > f-i + [i~ events according to the SM 
branching fractions, including cross-feed between the two 
modes. 

The contribution of the exclusive background compo- 
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TABLE I. Expected and observed limits on the B° — > 
fi + n~ branching fractions for the 2012 and for the combined 
2011+2012 datasets. 



Dataset 


Limit at 


90 % CL 


95 % CL 


2012 


Exp. bkg+SM 
Exp. bkg 
Observed 


8.5 x 10~ 1() 

7.6 x 10" 10 
10.5 x 10~ 10 


10.5 x 10~ 10 
9.6 x 10~ 10 
12.5 x 10" 10 


2011+2012 


Exp. bkg+SM 
Exp. bkg 
Observed 


5.8 x 10~ 1U 
5.0 x 10" 10 
8.0 x 10~ 10 


7.1 x 10~ 10 
6.0 x 10~ 10 
9.4 x 10~ 10 



nents is also evaluated for the 2011 dataset, modifying 
the number of expected combinatorial background in the 

fi + fi~ branch- 



signal regions. The results for the B^ 
ing fractions have been updated accordingly. We ob- 
tain B(B° -+ m + M~) < 5.1 x 1CT 9 and B(B° -+ fi + fj,~) 

< 13 x 10~ 10 at 95 % CL to be compared to the published 
limits B{B° S -+ n+^r) < 4.5 x 1(T 9 and B{B° -+ p+p,-) 

< 10.3 x 10- 10 at 95% CL 0, respectively. The (l-CL b ) 
p- value for B 9 — » changes from 18% to 11 % and 
the B 9 — > fi~ branching fraction increases by ~ 0.3 a 
from (0.8±i:|) x 10" 9 to (1.4±^) x KT 9 . This shift 
is compatible with the systematic uncertainty previously 
assigned to the background shape [§]. The values of the 
B 9 -+ branching fraction obtained with the 2011 
and 2012 datasets are compatible within 1.5 a. 

The 2011 and 2012 results are combined by computing 
the CL S and performing the maximum-likelihood fit si- 
multaneously to the eight and seven BDT bins of the 2011 
and 2012 datasets, respectively. The parameters that 
are considered 100 % correlated between the two datasets 
are f s /f d , B(B+ -> J/if>K+) and B{B° -> K + iz~), the 
transition point of the Crystal Ball function describing 
the signal mass lineshape, the mass distribution of the 
B 9 ^ -+ h + h'~ background, the BDT and mass distri- 
butions of the B° -+ Tr-fi+is^ and B°W -> 7r°W/i + /i- 
backgrounds and the SM predictions of the -B 9 — > 
and B° — > branching fractions. The distribution of 

the expected and observed events in bins of BDT in the 
signal regions obtained from the simultaneous analysis of 
the 2011 and 2012 datasets, are summarized in Table [TT] 

The expected and observed upper limits for the B° — > 
fi + n~ channel obtained from the combined 2011+2012 
datasets are summarized in Table |T] and the expected 
and observed CL S values as a function of the branching 
fraction are shown in Fig. [T] The observed CLb value 
at CL s+ b = 0.5 is 89%. The probability that back- 
ground processes can produce the observed number of 
B 9 -+ ^t + /i~ candidates or more is 5 x 10~ 4 and corre- 
sponds to a statistical significance of 3.5 a. The value of 
the B 9 — » branching fraction obtained from the fit 
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FIG. 1. CL S as a function of the assumed B° — > fi + fi~ branch- 
ing fraction for the combined 2011+2012 dataset. The dashed 
gray curve is the median of the expected CL S distribution if 
background and SM signal were observed. The shaded yellow 
area covers, for each branching fraction value, 34% of the ex- 
pected CL S distribution on each side of its median. The solid 
red curve is the observed CL S . 
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FIG. 2. Invariant mass distribution of the selected B° — >• 
fi + [i~ candidates (black dots) with BDT > 0.7 in the com- 
bined 2011+2012 dataset. The result of the fit is over- 
laid (blue solid line) and the different components detailed: 
B° — > /i + Ai~ (red long dashed), B° — > fi~ (green medium 



fJ- A 4 
dashed) , B 



(s) 



h + h'~ (pink dotted), B° -> n~ 



(black short dashed) and B 0(+) -> 7r 0(+ VV~ (light blue 
dot dashed), and the combinatorial background (blue medium 
dashed). 



is 



B{B\ 



>o 



' M + M ) : 



: (3.2ii; 4 (stat)iEJ-j(syst)) x 10~ 9 



and is in agreement with the SM expectation. The in- 
variant mass distribution of the B9s — > candidates 
with BDT > 0.7 is shown in Fig. |2J 



The true value of the B 9 — > branching fraction is 

contained in the interval [1.3, 5.8] x 10 _9 ([1.1, 6.4] x 10~ 9 ) 
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TABLE II. Expected combinatorial background, B? s s — > h + h! peaking background, cross-feed, and signal events assuming the 
SM prediction, together with the number of observed candidates in the _B° — > and B° — > mass signal regions, in 

bins of BDT for the 2011 (top) and for the 2012 (bottom) data samples. The quoted errors include statistical and systematic 
uncertainties. 

Mode BDT bin 0.0 - 0.25 0.25 - 0.4 0.4 - 0.5 0.5 - 0.6 0.6 - 0.7 0.7 - 0.8 0.8 - 0.9 0.9 - 1.0 

B° -+ Exp. comb, bkg 1880i| 55.5l 3 2 ;° 12.1±i;| 4.16t° : ?g ^ + -°ofi 0-77t° f s 0.47l ;f 6 0.24l°,;« 

(2011) Exp. peak, bkg 0.13l° °? 0.071 ,; , 2 . O.OStg;™ 0.05±g;g? 0.051 ,; , 2 0.05l° ; 02 0.05l ; 02 0.051°:°? 

Exp. signal 2.70t°; 8 80 ^°oZ l-03tg;f 7 0.92t°;l§ 1.06*>% l.W±° .{l 1.261°:?? ^±°of 5 

Observed 1818 39 12 6 1 2 1 1 

B° ^ mV 1 Exp. comb, bkg 1995+^ 59.2 ' r! 12^ 4.441°$ L67l^ 0.75±g : |g 0.44±°;§* 022+^ 

(2011) Exp. peak, bkg 0.781 ,; 3 * OAO^it 0.31±g".J| 0.281°,:™ 0.3ll° ; ° 8 0.30l ; ° 7 0.3ll°,; ° 0.30l° ; 08 

Exp. cross-feed 0.431°,;" 0.21±g;g* 0.161 ,; , 3 0.15±g;gl 0.17l°,;°| 0.17±g;gi 0.20l°,;°| 0.21±g;g| 

Exp. signal 0.331°,;!° 0.16l°;°I 0.131°;°^ 0.11±g : ™ 0.131 ,; , 2 0.13l ;° 2 0.15l°;° 2 2 0.16±g : g| 
Observed 1904 50 20 5 2 1 4 1 

Mode BDT bin 0.0 - 0.25 0.25 - 0.4 0.4 - 0.5 0.5 - 0.6 0.6 - 0.7 0.7 - 0.8 0.8-1.0 



B° -> 


Exp. comb, bkg 


2345l« 


56.7l 3 ;° 


13.llii 


4.421°;!! 


2.iol° ;^ 


S5+°- 42 

U.OJ_ 22 


0.39l°|? 


(2012) 


Exp. peak, bkg 


0.2501°;°? 


0-151°;°! 


0.08l°;° 3 2 


o.o8t°:£ 


0.07l°;° 2 2 


06+ 02 

U.UD_ 02 


in+ 03 




Exp. signal 


3.69l°;g 


2 14+ - 37 


i-2o±g : ?J 


1 16+ - 18 

1 - 1D -0.16 


17 +0.18 
J^- 1 ' -0.16 


I 1CC+0.19 
1 - lo -0.17 


9 i q+0.33 
^• lo -0.29 




Observed 


2274 


65 


19 


5 


3 


1 


3 



B° Exp. comb, bkg 249llg 59.51 3 ; 3 13.9±i;g 4.74lJ;°° 2.10l°;£ 0.55l°;|° 0.29l°; 34 

(2012) Exp. peak, bkg 1.49l°; 5 3 ° 0.86l° ;g 0.481°,;!°. 0.441°,;!? 0.42l°;^ 0.37l°; 3 0.62l°^ 

Exp. cross-feed 0.63l°; ° 0.36l° °? 0.20±g : g* 0.20l°;if 3 0.20±g : g| 0.20l°,;°| 0.361°;°° 

Exp. signal 0.441°;°° 0.26±g : g* 0.14±g-™ 0-1^ 0.14l°;° 2 0.14l°;° 2 0.26±g : g| 

Observed 2433 59 19 3 2 2 2 



at 90% CL (95% CL), where the lower and upper limit 
are the branching fractions evaluated at CL s+b = 0.95 
(CL s+b = 0.975) and CL s+b = 0.05 (CL s+b = 0.025), 
respectively. These results are in good agreement with 
the lower and upper limits derived from integrating the 
profile likelihood obtained from the unbinned fit. 

In summary, a search for the rare decays B® — > 
and B° — > /i + yU~ is performed using l.Ofb -1 and l.lfb" 1 
of pp collision data collected at y/s = 7 TeV and y/s = 8 
TeV, respectively. The data in the B° search window 
are consistent with the background expectation and the 
world's best upper limit of B(B° -> < 9.4 x 10~ 10 

at 95 % CL is obtained. The data in the B® search 
window show an excess of events with respect to the 
background-only prediction with a statistical significance 
of 3.5cr. A fit to the data leads to B(B® — > /i + /i~) 
= (3.2 j^ij) x 10 -9 which is in agreement with the SM 
prediction. This is the first evidence for the decay 
B° s ^»+»-. 
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